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METHOD AND APPARATUS FOR EVALUATING ERROR CONTROL 
PARAMETERS OF SELF-SIMILAR CONSTELLATIONS 

The present application is related to application serial 
number 09/741,636 filed on December 20, 2000, entitled "Performance Evaluation of 
Multi-Carrier Channels, " by inventors Vlad Mitlin and Richard G.C. Williams (Atty. 
Doc. No.: 3Com 72-l(3292TDCUSP) ) which is incorporated herein by reference in 
its entirety. 

FIELD OF THE INVENTION 

The present invention relates to multi-carrier data communications; 
and, more particularly, to evaluating the error control parameters of self-similar 
constellations. 

BACKGROUND OF THE DISCLOSURE 

Demand for digital subscriber line (DSL) service across existing 
twisted pair copper wires between a central office and a remote location is increasing. 
Typically, DSL services operate in accordance with DSL standards recommended by 
the Telecommunication Standardization Sector of the International 
Telecommunication Union (ITU). A family of DSL Recommendations from the ITU 
includes: G.992.1, G.992.2, G.991.1, G.996.1, G.994.1, G.997.1 and G.995.1. 
Recommendation G.995.1 provides an overview of these standards. 
Recommendations G.991.1, G.992.1, G.992.2 have developed techniques for 
transmitting a range of bit rates over the copper wires of the local network including 
high bit rates at relatively short distances, and lower bit rates at longer distances. In 

particular, the G.992.1 and G.992.2 recommendations are based on asymmetric 
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digital subscriber line technology that has different data rates in each direction of 
transmission. The G. 992.1 recommendation is referred to as G.dmt and uses a splitter 
to filter the voicegrade signals at the remote location. The G.992.2 recommendation 
is referred to as G.lite and does not use a splitter. Recommendations G.994.1, 
G.996.1 and G.997.1 support the G.992.1 and G.992.2 recommendations by providing 
common handshake, management and testing procedures. These standards allow 
substantial flexibility in implementation. 

DSL services typically use a discrete multi-tone (DMT) signal to 
transmit data. A DMT signal has multiple sub-channels, each of which is assigned a 
frequency, also referred to as a carrier frequency or a tone, belonging to a discrete 
frequency band. Because individual sub-channels operate at different frequencies, the 
sub-channels may have different operating characteristics. For instance, more power 
may be used at higher frequencies. In addition, different numbers of bits may be 
loaded on different sub-channels in accordance with their capacity, which depends on 
frequency, power, signal-to-noise ratio and transmission line characteristics. 
Sub-channels that do not meet or exceed a minimum signal-to-noise ratio are not 
used. When initiating a DSL communication session, in the DSL modem, an 
initialization procedure at the receiver determines a number of bits to be grouped into 
a symbol for each sub-channel, that is, a number of bits per sub-channel, and 
exchanges that information with the transmitting DSL modem. 

Quadrature amplitude modulation (QAM) is a technique to encode 
multiple bits into a QAM symbol. Each QAM symbol represents a distinct 
combination of bit values using a distinct combination of amplitude and phase of the 
carrier waveform. Each QAM symbol is represented by a QAM waveform. 
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Referring to Fig. 1, a signal space diagram depicts a constellation 14 of 
a group of distinct QAM symbols 18 that represents combinations of a group of bits. 
Fig. 1 depicts a 16-point G.992.2 constellation in which the number of bits b is equal 
to 4. In QAM, the amplitudes of two quadrature carriers are modulated and the 
5 carriers are combined. The x-axis 15 represents the amplitude of a first carrier, and 

the y-axis 16 represents the amplitude of a second carrier that is shifted in phase by 
90° with respect to the first carrier. For example, the first carrier is a sine wave, while 
the second carrier is a cosine wave. Each point 19 represents a distinct combination 
of the modulated carriers and thus a distinct QAM symbol. 

10 

A constellation encoder encodes groups of bits into QAM symbols. 
For example, for QAM symbols that represent four bits, the constellation will have 
sixteen distinct QAM symbols and map each of the sixteen possible combinations of 
the four bits to a distinct one of the QAM symbols. 

15 

In DMT systems, the digital information is transformed by a modem 
into an analog form that is essentially a sequence of DMT symbol waveforms. Each 
DMT symbol bears information in an array of zeroes and ones, which has several 
bi-sized sub-arrays. Each sub-array corresponds to a QAM waveform representing a 
20 2 bl -point constellation. In other words, b ( represents a number of bits per 

sub-channel z. A DMT symbol waveform is the superposition of these QAM 
waveforms. The channel itself is characterized by a signal-to-noise ratio y , where y 
represents the signal-to-noise ratio (SNR) at the i th carrier frequency. 

25 In DMT systems, each sub-channel has a constellation encoder. 

Typically, in DMT systems, equal error protection is applied, and the number of bits 
per sub-channel b t is determined as follows: for each sub-channel, the bit error rate 
should not exceed a target bit error rate prior to decoding and retransmission p b . One 
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conventional QAM error determination procedure as described by John G. Proakis, in 
Digital Communications, (Proakis) 1995, on p. 280 yields the QAM symbol error rate 
rather than the target bit error rate. Evaluating the equivalent bit error rate is known 
to be a complicated problem (See Proakis, p. 441). Therefore, it is typically assumed 
5 that the QAM symbol error rate is approximately equal to the bit error rate. However, 
this assumption is not precise and, for large QAM constellations, may misevaluate the 
bit error rate by a factor of -10. 

The assumption of an equally probable constellation decoding error 
0 yields an average fraction of erroneous bits that approaches Vi at large values of b. In 

other words, approximately half of the bits will be in error when a QAM symbol error 
occurs. This approach was earlier used in Proakis, p. 262 for 2*-ary orthogonal 
signals. However, this approach is not precise either. Because this assumption is 
used when determining the number of bits b t per sub-channel, an improved method 
5 and apparatus are needed to select the number of bits per sub-channel. Furthermore, 

this assumption does not accommodate for fluctuations in the bit error rate. The 
method and apparatus should also accommodate for fluctuations in the bit error rate. 

Reed-Solomon encoding is a method of forward error correction used 
3 in DSL communications to detect and correct transmission errors, effectively 

increasing the signal-to-noise ratio of the communications channel. By encoding 
information, errors may be reduced without decreasing the data rate. In 
Reed-Solomon encoding, redundant symbols are added to information symbols to 
allow errors to be detected and corrected. As the number of redundant symbols 
> increases, a greater level of noise may be tolerated. Among the transmission 

parameters to be selected in DSL communications are forward error correction (FEC) 
parameters for Reed-Solomon encoding. The FEC parameters determine the amount 
of information data in an information field, and a number of redundancy symbols that 
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are associated with the information field of an information frame. Reed-Solomon 
encoding parameters for forward error correction are selected while executing the 
initialization procedure. During initialization, the channel is analyzed and FEC 
parameters are determined and exchanged. The ITU recommendations set a bit error 

5 rate (BER) standard of 10" 7 . The Reed-Solomon encoding parameters for forward 

error correction are selected based on, at least in part, the bit error rate standard and 
an average number of erroneous bits per QAM symbol. However, there is a need for 
selecting FEC parameters based on a more precise determination of the average 
number of erroneous bits per QAM symbol. In addition, there is a need to select FEC 

10 parameters based on a target bit error rate that accommodates for fluctuations in the 
average number of erroneous bits per QAM symbol. 



SUMMARY OF THE INVENTION 



15 A method, and apparatus implementing the method, selects a bit load 

for a channel. The channel may be a single carrier or multi-carrier channel. In one 
aspect of the invention, the channel has a single carrier. Bits are encoded based on a 
constellation of points. The constellation has a self-similarity property. The bit load 
for the channel is selected based on the self-similarity property of the constellation. 

20 

In another aspect of the invention, a multi-carrier system encodes bits 
based on a constellation of points. Each point represents a tuple of bits. The 
constellation has a self-similarity property. The bit load per sub-channel is selected 
based on the self-similarity property of the constellation. 

25 

In yet another aspect of the invention, the constellation is square. 
Alternately, the constellation is non-square. In another aspect of the invention, the 



probabilistic properties of non-square constellations can be described by the same 
equations that are used for even-bit or square constellations. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 is a signal space diagram of an exemplary QAM constellation 
with a number of bits b equal to 4. 

Fig. 2 is a block diagram of a DSL modem at a central office, 
transmission medium, and DSL modem at a remote terminal. 

Fig. 3 is a block diagram of a transmitter of the DSL modem of Fig. 2. 

Fig. 4 is a block diagram of a receiver of the DSL modem of Fig. 2. 

Fig. 5 is a signal space diagram that depicts an exemplary constellation 
with a number of bits b equal to 2. 

Fig. 6 depicts the expansion of a constellation into the next larger 
square constellation. 

Fig. 7 depicts a construction used to prove proposition 2. 

Fig. 8 depicts a geometrical representation of error control properties 
of QAM constellation encoders using a quarter of a 256-point QAM constellation. 



Fig. 9 depicts a geometrical representation of error control properties 
of an odd-bit QAM constellation encoder, that is, a quarter of 32-point QAM 
constellation. 
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Fig. 10 is a general flowchart of a modem initialization module that 
determines a bit load based on a self-similarity property of a constellation in 
accordance with an embodiment of the present invention. 

5 Fig. 1 1 is a general flowchart of an alternate embodiment of a modem 

initialization module which determines a bit load based on the self-similarity property 
of a constellation and also accommodates for fluctuations in the number of bit errors 
per QAM symbol in accordance with the present invention. 

10 Fig. 12 is a more detailed flowchart of a modem and a channel 

initialization module which determine a bit load based on the self-similarity property 
of a constellation and accommodates for fluctuations in a number of bit errors per 
QAM symbol in accordance with an embodiment of the present invention. 

15 Fig. 13 is a more detailed flowchart of a channel initialization module 

which determines a bit load based on the self-similarity property of a constellation 
and accommodates for fluctuations in a number of bit errors per QAM symbol in 
accordance with an alternate embodiment of the present invention. 

20 To facilitate understanding, identical reference numerals have been 

used, where possible, to designate identical elements that are common to some of the 
figures. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

25 



A technique for determining the bit load of a channel is based on the 
error control properties of QAM constellation encoders. In particular, this technique 
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is applied to determining a bit load for a DMT or, which is sometimes referred to as 
an orthogonal frequency division multiplexed (OFDM), channel. 

Often, it is assumed that an error in constellation decoding yields any 
5 constellation point with the same probability. This assumption yields a value of l A for 
co, an average fraction of erroneous bits per an erroneous QAM symbol. However, in 
a bit-error-rate-controlled transmission environment this assumption is not true. It 
has been suggested that the main contribution to the error-control parameters of QAM 
constellation encoders is made by the nearest-neighbor errors. However, the 
10 expressions for the error control parameters of QAM constellation encoders may be 
cumbersome because direct calculations on QAM lattices are involved, and may be 
difficult to use directly in performance evaluation of multi-carrier channels. 

The present invention describes a technique to estimate the 
15 dependence of error control parameters of QAM constellation encoders with respect 
to the bit size of the QAM symbol This technique uses a self-similarity property of 
QAM constellations, that is, a distribution of a parameter in a sub-block of a 
constellation looks similar to the whole constellation, but on a smaller scale. In 
particular, in the present invention, the self-similarity property is the Hamming 
20 distances of the points in the constellation. Using the self-similarity property, the 

probability of having m bit errors in an erroneous b-bit QAM symbol is determined. 
Based on this result, an average fraction of erroneous bits o in an erroneous QAM 
symbol is determined and used to determine the bit load of a channel. When the 
number of bits of a QAM symbol b is large, the average fraction of erroneous bits co 
25 is approximately inversely proportional to b/2. A parameter called the mean square 
deviation measures fluctuations in bit error rate (BER) that are inherent to QAM 
constellation encoders. The dependence of the mean square deviation with respect to 
the size of the QAM symbol b in an erroneous QAM symbol is determined. The 
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method is developed for square constellations. In an alternate embodiment, 
non-square constellations are represented as a set of adjacent square constellations of 
different sizes, therefore reducing the determination of their error-control parameters 
to that of square constellations. When the constellation size is large, the probabilistic 
5 properties of non-square constellations can be described by the same equations that 
are used for square constellations. This approximation can be used in error control 
determination for single and/or multicarrier QAM channels where both square and 
non-square constellation encoders are applied, for example, G.992.2 (G.lite) and 
G.992.2 (G.dmt) asymmetric digital subscriber line (ADSL) standards. 

10 

1. An Exemplary DMT Communications System 

Referring to Fig. 2, a block diagram of an exemplary discrete 
multi-tone (DMT) communications system or channel 20 is shown. In one 

15 embodiment, the DMT communications system 20 is a DSL system. Alternately, the 

DMT communications system 20 is an asymmetric DSL (ADSL) system. In another 
alternate embodiment, the communications system 20 has a single QAM channel. 
The DMT communications system 20 connects to a core network 22 at a central 
office 24. The core network 22 may include functions such as a concentrator, a 

20 switch and interfaces to broadband and narrowband networks. The DMT 

communications system 20 has a modem 26 at the central office 24 and a modem 28 
at a remote terminal 30 that are interconnected by a transmission medium 32. The 
modems 26, 28 implement a multicarrier transmission method in accordance with the 
present invention. The modem 26 at the central office 24 includes a processor 

25 (CPU) 34, a transmitter/receiver (Tx/Rx) 36, and a memory 38, which are 

interconnected by a bus 39. The memory 38 stores a modem driver 40 that is 
executed by the processor 34, and tables 42 in accordance with various embodiments 
of the present invention. At least one table 42 stores values representing a number of 
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bits per sub-channel 46 for each of the sub-channels. One or more of the tables 42 
store FEC parameters. In one embodiment, the modem driver 40 includes a modem 
initialization module 48 that generates the various tables of the present invention; and 
a channel initialization module 43 that accesses the tables to select a bit load. In an 

5 alternate embodiment, the channel initialization module 43 is not executed once, but 

may be executed periodically to monitor the data, update the tables 42, and adjust the 
FEC parameters, if needed, and also to adjust the number of bits per sub-channel, if 
needed, during system operation. In another alternate embodiment, the channel 
initialization module 43 generates the number of bits per sub-channel and FEC 

10 parameters based on a measured effective signal-to-noise ratio, rather then the modem 

initialization module 48. The memory 38 may be implemented using RAM, ROM, or 
a disk drive, or a combination thereof. In an alternate embodiment, the modem driver 
procedure 40, and/or portions thereof, may be implemented in hardware using digital 
logic, for example, using an application-specific integrated circuit (ASIC). 

15 

The modem 28 at the remote terminal 30 includes the same 
components as the modem 26 at the central office 24, and also includes an external 
interface 50 to connect to a computer system. The external interface 50 may be a 
PCI, internet, or USB interface. In an alternate embodiment, the modem 38 is 
20 internal to the computer system. 

The communications channel uses a discrete multitone (DMT) signal 
to carry data downstream and upstream on the transmission medium 32. The 
downstream direction is from the central office 24 to the remote terminal 30; and the 
25 upstream direction is from the remote terminal 30 to the central office 24, Typically, 

higher data rates are supported for transmitting data in the downstream direction than 
the upstream direction. 
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In another embodiment, the method is applied to DSL systems 
compliant with the ITU Recommendation G.992.2 standard. The G.992.2 standard 
specifies allowable sets of Reed-Solomon code parameters that are applied equally 
across all sub-channels. The optimal coding parameters are found depending on, at 
least in part, the signal-to-noise ratio (SNR) distribution over the channel's frequency 
band and on the number of sub-channels used. Alternately, the present invention is 
applied to a DSL system compliant with the ITU Recommendation G.992.1 standard. 
In an alternate embodiment, the present invention can be applied to any multi-carrier 
system. 

In another embodiment, the transmission medium is twisted pair. 
Alternately, the transmission medium may be a coaxial cable. In another alternate 
embodiment, the transmission medium is "wireless," that is, the transmission medium 
has radio transmitters and receivers to transmit and receive the DMT signal, 
respectively. Alternately, the transmission medium is implemented optically. In one 
embodiment, the transmission is performed using an infrared (IR) interface which 
includes an IR transmitter and an IR receiver. Alternately, the transmission medium 
is implemented using optical fiber with an optical transmitter and receiver. 

Fig. 3 is a block diagram of a transmitter 60 of the transmitter/receiver 36 
of Fig. 2. A scrambler and FEC block 62 scrambles an incoming stream of data bits 
using a predetermined algorithm, and applies FEC to the bits to produce FEC encoded 
bits on lead 66. A tone ordering block 68 assigns the FEC encoded bits to the usable 
sub-channels and outputs the assigned FEC encoded bits on leads 72. Constellation 
encoders 74, one for each sub-channel, convert the FEC encoded bits assigned to a 
sub-channel to a QAM symbol or constellation point which are output on leads 76. An 
inverse discrete Fourier transform (IDFT) block 78 modulates the output from the 
constellation encoders 74 to produce a digital modulated signal representing the 
superposition of the QAM symbols on lead 80. A parallel-to-serial output buffer 82 
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stores the digital modulated signal from the IDFT block 78. A digital-to-analog converter 
(DAC) 84 serially converts the stored digital modulated signal from lead 86 to an analog 
DMT waveform that is output on lead 88 for transmission. The functions of the 
blocks 62, 68, 74, 78, 82 and 84 may be implemented completely in hardware as a circuit, 
5 or alternately, using a combination of hardware and software. 

Fig. 4 is a block diagram of a receiver 90 of the transmitter/receiver 36 of 
Fig. 2. The receiver 90 performs the inverse of the functions of the transmitter 60 to 
produce a decoded bit stream from the analog DMT waveform. An analog-to-digital 

10 converter 92 receives the analog DMT waveform on lead 94 and converts it to a series of 
digital values to produce a digital DMT waveform on leads 96. A serial-to-parallel input 
buffer 98 stores the values of the digital DMT waveform and outputs the values, in 
parallel, on leads 102 to the discrete Fourier transform block 104. The discrete Fourier 
transform block 104 demodulates the values of the digital DMT waveform to produce a 

15 set of QAM symbols, one for each usable sub-channel, on leads 106. Constellation 

decoders 110 decode their respective QAM symbols and output respective sets of FEC 
encoded bits for the usable sub-channels on leads 1 12. A tone reordering block 1 14 
places the FEC encoded bits in their proper order and outputs the ordered encoded bits on 
leads 116. The descrambler and FEC decoder 120 corrects for errors using forward error 

20 correction and descrambles the bits to reproduce the data stream sent by the transmitter 

on lead 122. 

2, QAM constellation encoders 

25 In multicarrier systems, the digital information is transformed by the 

modem into an analog form that is essentially a sequence of DMT symbol 
waveforms. Each DMT symbol bears information in an array of zeroes and ones, 
which has several bi-sized sub-arrays. Each sub-array corresponds to a QAM 
waveform representing a 2 bl -point constellation. In other words, b t represents the 
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number of bits per sub-channel i. A DMT symbol waveform is the superposition of 
these QAM waveforms. The channel itself is characterized by a signal-to-noise 
ratio Yu where y t represents the signal-to-noise ratio (SNR) at the i th carrier frequency. 



5 An average fraction of erroneous bits in an erroneous 6,-sized QAM 

symbol is represented by co(y u h). If pqam(Yu h) is a symbol error rate, the bit error 
rate is equal to the product of the symbol error rate and the average fraction of 
erroneous bits in an erroneous 6 r sized QAM symbol, that is, Pqam(Yi > bi) ®( Yb bt). 
For a multicarrier system, for each sub-channel, the bit error rate (BER) is less than or 

10 equal to the target bit error rate pb in accordance with the following relationship: 

BER = Pqam (r, > b, M 7i , Z> f ) < />> , (2.1) 

where p b is the largest BER value allowed in the system past the constellation 
15 decoding stage. 

In one embodiment, in the absence of a forward error correction block 
in a transceiver, the G.lite and G.dmt standards specify that the value of the target bit 
error rate/?* is equal to 10~ 7 . To improve performance, the number of bits per 
20 sub-channel b t should maximize the left hand side of Equation (2.1). 

While an expression for the symbol error rate pqam is well known from 
communication theory: 



25 P QA m 



(2.2) 
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an appropriate expression for the average fraction of erroneous bits in an erroneous 
6 r sized QAM symbol mfyu bt) will be determined. The problem is formulated as 
follows: assuming that a 6-sized sub-channel is in error, determine the average 
fraction of bits in error in such a sub-channel, co(y,b). 

One way to approach this problem is to assume that if a QAM symbol 
is received in error, that symbol will result in obtaining any constellation point with 
the same probability. For a Z>-bit QAM constellation, this yields: 



16 II ryb-X 

a(p) = iy m ^fL = ^_ (2.3) 



Equation (2.3) follows from the fact that the probability of having mi's (m>0) in a 

h 

randomly chosen binary array of the length b is equal to /(2 - 1) . 

w 



Equation (2.3) shows that the assumption of an equally probable 
constellation decoding error yields an average fraction of erroneous bits in an 
erroneous symbol o)(b) that approaches l A at large values of &. In other words, 
approximately half of the bits will be in error when a QAM symbol error occurs. 
This approach was earlier used in Proakis, 1995, p. 262 for 2^-ary orthogonal signals. 

However, in practice, in a bit-error-rate-controlled transmission 
environment, the results of decoding an erroneous symbol or point in the constellation 
tend to group around the correct, that is, the transmitted, constellation point. Usually 
the results of decoding the erroneous symbol are the nearest neighbors of the correct 
point on the constellation. For QAM systems, if only the nearest neighbors of a point 
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in a constellation are considered, the problem of computing the average fraction of 
erroneous bits per symbol co(b) becomes geometrical. 

Referring back to Fig. 1, in the exemplary 16-point G.992.2 
constellation 14, a coordination number represents a number of nearest neighbors to a 
point. For example, point 6 has two nearest neighbors, points 4 and 7. More 
generally, for an arbitrary 2 6 -point square constellation, the coordination number is 
equal to 4 for each internal point, 3 for each external non-corner point, and 2 for each 
external corner point of a constellation. 

Let ai be the binary representation of the label for the z th point of the 
constellation. Let %i ^ e the coordination number of this z th point, that is, the total 
number of nearest neighbors to the i ih point in the constellation. Then the average 
fraction of erroneous bits per symbol co(b) is determined in accordance with the 
following equation: 




In Equation (2.4), d H (; •) is the Hamming distance between two tuples, and the 
innermost summation in Equation (2.4) is performed over the nearest neighbors of the 
tuple or point i. It has been found that l/co(b) is an increasing function of the number 
of bits per sub-channel b. At large values of b, l/co(b) increases linearly with a 
proportionality coefficient of the order of V%. The following relationship for l/o)(b) 
has been suggested: 



l/o)(b) = (2b + 3)/4 
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Equation (2.5) yields the asymptotic behavior discussed above. 

Since Equation (2.5) is one possible approximation of Equation (2.4), 
it is important to determine the average fraction of erroneous bits per symbol co(b ) 
based on Equation (2.4). However, solving Equation (2.4) is cumbersome and, for 
large values of b, is difficult to use in direct calculations of error control parameters. 
A more general problem is to develop a method of determining other error control 
parameters of QAM constellation encoders, that is, determine the probability of 
having m bit errors in an erroneous 6-bit QAM symbol. Therefore, the geometrical 
structure of QAM constellations that directly affects the error control properties of the 
constellation encoders will be discussed. Initially, square QAM constellations, in 
which the number of bits b is even, will be discussed. 

For a given QAM channel, or alternately, in DMT communications, 
for a given QAM sub-channel, the constellation encoder selects an odd-integer 
point (X, Y) from the square-grid constellation based on the binary 6-bit array, or 
tuple, (v h „ h v h -2, ..,vi,vo}. For convenience of description, these b bits are identified 
with an integer label whose binary representation is {v b . 1} v b . 2f ..,v;,v^. 

Referring to Fig. 5, for example, for b=2, the four constellation points 
are labeled 0, l,2and3. The integer values X and Y of the constellation point are 
determined from the b bits of the input tuple {v b . h v b . 2f ..,v 2 ,v 0 } as follows. Xand Y 
are the odd integers with 2's complement binary representations {v b . I} v b . 3 , .,,v L 1} 
and {v b - 2 ,v h - 4j ..,v 0 , 1}, respectively. For example, when b=2, and the input bit 
stream {v h v 0 } is equal to {1, 0}, X={1,1} andY={0,l}. 

Referring to Fig. 6, a constellation can be expanded into the next larger 
square constellation. For example, the 4-bit constellation of Fig. 1 can be obtained 
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from a 2-bit constellation by replacing each label n by a 2 x 2 block of labels. The 
same procedure can be used to construct larger even-bit constellations recursively. 
The constellations obtained for even values of b are square in shape. 

5 3. Self-similarity of error control parameters of QAM constellation encoders 

A recursive procedure of constructing larger QAM constellations from 
smaller ones forms the basis for analyzing their error control properties. The 
following model will be used: assuming that the digital-to-analog conversion yields 
10 constellation point / at the transmitter, an error of the constellation decoder at the 
receiver will yield one of the nearest neighbors of i. Also, assume that a decoding 
error results, with equal probability, in any of the nearest neighbors of the 
constellation point i transmitted. Therefore, this probability is equal to one over the 
coordination number of the constellation point L 

15 

The following proposition will be used. 

Proposition 1 . 

20 If % b, c, d, and n are integers such that c<2 n , d<2 n , and d H (..,..) is the 

Hamming distance between the binary representations of two integers, then 

d H (2 n a+c, 2 n b+d)= dn(a,b)+ d H (c,d) (3.1) 

25 Proof. 



Binary representations of 2 n a and 2 n b have zeroes in their n least 
significant bits. The total number of bits in binary representations of c and d does not 



-18- 

exceed n. Therefore, Equation (3.1) follows directly from the definition of the 
Hamming distance, e.g., the total number of different bits in two binary arrays. 

The following proposition will also be used. 

Proposition 2 . 

Let a and b be the labels of 2 neighbor points with the Hamming 
distance of d h (a,b) in a 2n-bit QAM constellation. Let a (2n+2)-bit constellation be 
generated from a 2n-bit constellation according to the rule shown in Fig. 6. Then the 
Hamming distance between two neighbor points of two different 2x2 blocks 
generated out of a and b is equal to d h (a,b)+l . 

Proof 

Referring to Fig. 7, any two neighbor points of two different 2x2 
blocks 140 and 142 generated out of a and b have the same parity, that is, either both 
are even, or both are odd. As Fig. 7 depicts the case of b 142 being on the right of 
a 140, it can be seen that this is true for b being on the left, above, and below a. 
Applying Equation (3.1) proves the proposition. 

Referring to Fig. 8, propositions 1 and 2 allow the error control 
properties of QAM constellation encoders to be described geometrically. Fig. 8 
depicts a quarter of an 8-bit a constellation that has 256 constellation points. Small 
squares 152 represent the constellation points. Hamming distances of 1, 2, 3 and 4 
between the adjacent constellation points are shown by arrows 154, 156, 158 and 160, 
respectively. Shaded crosses, indicated by the arrows, represent the boundaries 
separating any two points with the Hamming distance of n. For instance, if a 
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decoding error yields a point lying across the 4-boundary, the incorrect decoding 
yields 4 bit errors; across the 3-boundary yields 3 bit errors; across the 2-boundary 
yields 2 bit errors; and across the 1 -boundary yields 1 bit error. The constellation of 
Fig. 8 has one 4-cross having a linear size of 16, four 3-crosses having a linear size of 
8, 4 2-crosses having a linear size of 4, and 4 1 -crosses having a linear size of 2. 

Another way to consider the error control properties of QAM 
constellations is to examine the distribution of Hamming distances along a horizontal 
or vertical line across a constellation. In the 2 2 ~point constellation, the distribution of 
Hamming distances is: 

1. 

In the 2^-point constellation the distribution of Hamming distances is: 
1 2 1. 

In the 2^-point constellation the distribution of Hamming distances is: 

12 13 12 1. 

Generally, the distribution in a larger 2 b+2 -point square constellation can be generated 
from the distribution in a smaller 2*-point constellation by increasing the values in the 
latter distribution by one and then inserting Ts at the beginning, between each two 
neighbors, and at the end. 

The distribution of Hamming distances repeats itself inside smaller 
blocks of the constellation. There is a central cross in each square block, and the 
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constellation neighbors separated by the cross have the maximal Hamming distance 
in the block. Each of the four sub-blocks of the square block with the n-error central 
cross repeats the pattern having an (n-l)-error central cross. Such a geometrical 
structure is called self-similar. 

5 

Examples of self-similar structures are well-known in modern fractal 
theory, the hydrodynamic analysis of turbulence, the theory of critical phenomena, 
and the mechanics of composite materials. Self-similarity can be considered as a kind 
of advanced symmetry. In the present invention, the self-similar properties of QAM 
10 constellations are used to determine the error control parameters for constellation 
decoders. 

4. Analysis of the error probability distribution of QAM constellation encoders 

15 Referring to Fig. 8, consider a square 6-bit QAM constellation 

represented by a 2 N x 2 N square; thus the number of bits b is equal to 2N, b=2N. The 
coordination number %i of the i th point is equal to the total number of nearest 
neighbors to the i th point in the constellation. Using the nearest-neighbor-error 
supposition above, the probability of having k bit errors in an erroneously decoded 

20 QAM symbol or tuple is determined as follows: 

^ i=l d H (a l> a J )=k Xi 

An (i ,y)-link represents an edge starting at the constellation point i and 
25 ending in a nearest neighbor point j. The weight of a (z j)-link is equal to one over the 
coordination number (Vzd- The maximum number of bit errors that can be made in 
an erroneously decoded tuple, under the nearest-neighbor-error supposition, is equal 
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to N. The links intersecting the central cross of the f 1 x f 1 square contribute to the 
probability of having N bit errors in an erroneously decoded tuple p(N 9 b) in 
accordance with the following relationship: 



iN,b): 



>2iV 



1 (\ \\ 

1-4-2^+8 i-i 1 
4 13 4 



(4.2) 



Equation (4.2) accounts for the fact that 8 links corresponding to 4 end points of the 
cross have a weight of 1/3, not 1/4. 

10 Similarly, the probability of having N-l bit errors in an erroneously 

decoded tuple p(N-l,b) is contributed by all the links intersecting the central crosses 
of the four adjacent x 2 1 * 1 ' 1 squares forming the 6-bit constellation in accordance 
with the following relationship: 



15 



P 



(N-l,b)- 



>2JV 



1 ( 

4---4-2"" 1 +2-8 



3 4. 



(4.3) 



Equation (4.3) accounts for the fact that 16 links corresponding to 8 end points of the 
4 crosses have the weight of 1/3, not 1/4. 
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Generally, the probability of having N-i bit errors in an erroneously 
decoded tuple p(N-i,b) is contributed by all the links intersecting the central crosses of 
the 4 adjacent squares with the sizes f~ l x f' 1 forming the 6-bit constellation in 
accordance with the following relationship: 



25 



+ 8 



3~4 



0<i<N-2 



(4.4) 
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Equation (4.4) describes the probabilities of having k=N-i bit errors in 
decoding a QAM symbol, for k>L At k=l, the probability is determined in 
accordance with the following relationship: 



N-2 



P(W + T,p{N-i,b) = p{\,b) + -L g 



i=0 



2 N+t +S 



V3 4y 



2' 



-1 (4-5) 



Performing the summation in Equation (4.5) yields: 



\2N 



\2N-1 .4 ^ N _y 2 



+ --2 N ~ X 
3 



+ — 
3 



(4-6) 



Finally, noting that in the above equations i=b/2-k and N=b/2 yields: 



1 + - 



3-2 



6/2 



J_2 
+ 2* 3 



1 ^2_ 

2*/2 + 2* 



<y it , \<k<b/2 



(4.7) 



l,w/?ere/ = y] 



Equation (4.7) shows that at large values of b, the probability of having & bit errors in 
decoding a b-sized QAM symbol is determined in accordance with the following 
relationship: 



\+o 



sbll 



\<k<bl2. 



(4.8) 



and is asymptotically independent of b. In other words, for large constellations that 
have large values of b, the probability of having A: bit errors approaches l/2 k . 
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5. Average number of bit errors in an erroneous tuple 

Equation (4.7) allows the error control parameters of square QAM 
constellation encoders to be calculated. First, an average number of bit errors in an 
erroneous tuple <n e > is determined by applying Equation (4.7) to yield: 



b/2 



2 k if l 2 ^ 



£=1 



)b/2 n b 



(5.1) 



Since the sum on the right hand side of Equation (5.1) is equal to: 



V 2" 2*2"; 



(5.2) 



Equation (5.1) can be rewritten in the following form: 



J b 1 ^ ]_ 

2 hn 4 2 1 " 2 )3 



1 _2_ 

2 bl2 + 2 b 



= 2- 



J 2 b b _12-2" -(3b + 2)2 b/2 -2b-4 



(5.3) 



3-2 bn 3-2* 2-2* /2 3-2 



6-2 b 



At large values of b, the average number of bit errors in an erroneous tuple <n<> tends 
to 2. Therefore, setting co(b) to equal <n e >/b, the average number of bit errors per 
erroneous symbol co(b) is equal to: 



(b) 



_12-2*-(3ft + 2)2* /2 -26-4 
6b-2" 



(5.4) 
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Equation (5.4) shows that at large values of b, the value of l/co(b ) tends to b/2. This 
is what was observed when computing the average number of bit errors per 
symbol co(b) directly from Equation (2.4) and in numerical simulations of a QAM 
channel, and is contrary to the value of 2 discussed above. 

6. Mean square deviation of the number of bit errors in an erroneous tuple 

Similarly, higher moments of p(k,b) and the average number of bit 
errors in an erroneous tuple <n<> can be determined. In particular, the second 
moment of the average number of bit errors in an erroneous tuple <n e >will be 
determined, and an expression for the mean square deviation of the number of bit 
errors in an erroneous tuple, cr e9 will be derived. 

First, the second moment of p(k) is determined in accordance with the 
following relationship: 

Jfc=i V )h=\ £ ->\2 2 ) 

It is shown, below, that 

^k 2 n 2 +4n + 6 

LTk= 6 ^ • (6-2) 

Using Equation (6.2), Equation (6.1) can be rewritten as follows: 

/ 2\ f, 2 1 b 2 +Sb + 24] if 1 2} 
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Combining Equations (6.3) and (5.4) yields the final expression for <j e : 



V 



1+- 



3-2 



b/2 



b 2 +Sb + 24 



4-2 



b/2 



+ - 

3 



*b/2 



+ - 



12-2* -(3l> + 2)2* /2 -26-4 
6-2 b 



\ 2 



3b 2 +246 + 20 t (66 + 4)2* /2 -6 2 ~4 , f (36 + 2)2* /2 +26 + 4 



12-2 



6/2 



6-2* 



+ 



3-2' 



At large values of b, a e tends to the square root of two. 



(6.4) 



10 



Equation (6.4) yields a measure of the bit error rate fluctuations that 
are associated with using QAM constellation encoders. Therefore, a e determines the 
relative range of the BER fluctuations past the constellation decoding stage, i.e.: 



a(ber) 

BER 



(6.5) 



15 



Table 1, below, presents at least a portion of the error control 
parameters of square QAM constellation encoders for values of b varying between 2 
and 32. As shown in Table 1, l/o)(b) increases linearly at large values of b. Both 
<n e > and <r e 2 monotonically increase, approaching the value of 2 at large b. 
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Table 1. Error control parameters of QAM 6-bit constellation encoders. 



0 


co(b) 




<n e > 


2 

a e 


a e 


2 


5.000000e-00l 


2 


1 


0 


0 


A 

4 


3.229 167e-00l 


3.096774e+000 


1.291667e+000 


2.065972e-001 


4.545297e-001 


0 


f> r- y c\ AAA f\f\ -t 

2.569444e-001 


3.891892e+000 


1.541667e+000 


5.190972e-001 


7.204840e-001 


8 


2.l45l82e-00l 


4.661608e+000 


1.716146e+000 


8.543226e-00l 


9.242958e-001 


10 


l.829427e-00l 


5.466192e+000 


1.829427e+000 


l.l 6231 le+000 


1.078105e+000 


12 


l.583252e-00l 


6.3161 14e+000 


1.899902e+000 


L416250e+000 


1.190063e+000 


14 


1.387416e-001 


7.207642e+000 


1.942383e+000 


1.609229e+000 


1.268554e+000 


1 


1.229598e-001 


O 1 1^^7/1 1 ~ i AAA 

8. 13274 le+000 


1.967356e+000 


1.74697 le+000 


1.321730e+000 


18 


1.100970e-001 


9.082902e+000 


l.981745e+000 


1.840604e+000 


1.356689e+000 


20 


9.949509e-002 


L005075e+001 


l.989902e+000 


1.901852e+000 


1.379077e+000 


22 


9.065747e-002 


L103053e+001 


1.994464e+000 


L940705e+000 


1.393092e+000 


24 


8.320785e-002 


1.201810e+001 


1.996988e+000 


1.964748e+000 


1.401695e+000 


26 


7.686047e-002 


1.301059e+001 


l.998372e+000 


1.979325e+000 


1.406885e+000 


28 


7.139733e-002 


1.400613e+001 


L999l25e+000 


1.988016e+000 


1.409970e-K)00 


30 


6.665107e-002 


1.500351e+001 


1.999532e+000 


1.993123e+000 


1.411780e+000 


32 


6.249221e-002 


1.600199e+001 


l. 99975 le+000 


1.996089e+000 


1.412830e+000 



For example, fluctuations in the bit error rate may cause the bit error rate to exceed 
5 the target bit error rate. Therefore, an adjusted target bit error rate can be detenriined 
using the standard deviation of Table 1 for a specified constellation size. To 
accommodate for fluctuations, Equation 6.5 is rewritten as follows: 



10 



A(BER) (BER'-BERo) 

-^eT = (T °= ber 0 • (66) 



where BER 0 represents the adjusted target bit error rate, and BER' represents the 
unadjusted target bit error rate. For example, for a constellation with a number of 
bits b equal to 32, the mean square deviation o e is equal to 1.41283 and a target bit 
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error rate BER' of 10~ 7 , solving Equation (6.6) yields an adjusted bit error rate BER 0 
of 0.415 x 10" 7 . The adjusted target bit error rate BER 0 may then be used in 
subsequent calculations to determine FEC parameters. In this way, fluctuations in the 
bit error rate are accommodated, and the system error rate will meet the unadjusted 
target bit error rate. 

The derivation of Equation 6.2 will now be described. Consider the 
following function: 




(6.7) 



Rewriting Equation (6.7) yields: 



f(x) = x 




x 



2x 2 



(6.8) 




2"(x-2) 



x-2 



Differentiating Equation (6.7) twice, the following equality holds: 



/"(l) = S 2 +3E 1+ 2E 



(6.9) 



where 



A=0 



(6.10) 
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As 2 0 =2-2-" (6.11) 

=2-2-"(n + 2) (6.12) 



and 



ra) = 16 _fct^|+4)±2 (613) 



combining Equations (6.9) to (6.13) yields Equation (6.2). 
7. Non-square self-similar constellation encoders. 



The method of estimating error-control parameters developed above 
for square QAM constellation encoders is extended to non-square constellations 

15 which have an odd number of bits 6. A non-square constellation can be represented 

as a set of adjacent square sub-constellations of different sizes, therefore reducing the 
determination of its error-control parameters to the basic case considered in the 
present invention. Fig. 9 presents an upper left quarter 170 of a 32-point QAM 
constellation used in G.lite/G.dmt ADSL standards. Hamming distances of 1, 2, and 

20 3 between adjacent constellation points are shown by arrows. 

The entire constellation can be described as a 1 6-point QAM 
sub-constellation shown in Fig. 1, with 16 adjacent points over the perimeter of the 
1 6-point QAM constellation. The next odd-bit constellation, a 128-point QAM 
25 constellation, is generated from the 1 6-point constellation following the rule shown in 
Fig. 6. The rule of Fig. 6 may be repeatedly applied to generate larger constellations. 
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Therefore, the larger constellations have one central 2 2N -point QAM sub-constellation 
and 16 adjacent 2 2N ' 4 -point QAM sub-constellations. 

For a general, not necessarily square QAM constellation, the 
5 probability of having k bit errors in an erroneous QAM symbol is determined in 

accordance with the following relationship: 



p(k,b x = ^^j- p[k\^symbol transmitted e Qj & symbol received e Q ; ) 

7=1 2 



10 



15 



J 

+E 

j,m 



p{symbol transmitted e Q y & symbol received eQ m ) 
x j^Ajsymfo/ transmitted e Qj & symbol received e Q m ) 



(7.1) 



where 2 6j is the number of constellation points in each of /adjacent square 
sub-constellations Qj forming the non-square constellation considered, and 
sub-constellation Q is different from sub-constellation Q m \ and the following 
condition holds: 



7=1 ^ 



Notice that the Hamming distance between two neighbor points 
belonging to the same 2 bj -point square sub-constellation is not larger than bj/2<b/2. 
20 However, if two neighbor points of the constellation belong to different (adjacent) 

sub-constellations 5 the Hamming distance between them can exceed b/2. For 
instance, the Hamming distance between the points labeled by 9 and 19 in Fig. 9 is 
equal to 3. 



-30- 



Various non-square constellation encoders look like a set of adjacent 
square constellations of different sizes. As evaluating p(k,b 1 ,..,b J ) from 
Equation (7.1) in each specific case of a non-square constellation is quite 
cumbersome, the following general result can be derived: 

Proposition 3 . 

Let a non-square 2 b -point QAM constellation consist of J adjacent 
square sub-constellations. Let each 2 bj -point sub-constellation have been generated 
recursively according to the rule shown in Fig. 6. Denote 

b* =mmbj,j = l,.J . 

Then at large bj,j=l,..,J, and at 1 < A: < b, 12 ,p(k,b h ..,bj) becomes asymptotically 
close to p(k,b), the error probability distribution of a square constellation encoder 
(Equation (4.7)): 

p{k,b 1 ,..,b J )^p(k,b) , as bj -> coj = \,..J 

Proof. 

Consider different terms in Equation (7. 1). To determine the first 
conditional probability in the first sum on the right hand side of Equation (7.1), 
Equation (4.7) is used at 1 < k < b< 12 , or, alternately, Equation (4.8) is used at large 
values of b h for each sub-constellation. Generally, this conditional probability, as 
determined by Equation (4.8), is complemented by a correction term. This correction 
term accounts for the change in the coordination number of some external points of a 
square b r bi\ sub-constellation as a result of bringing it in contact with another 
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sub-constellation. For instance, a point with the value of 1 1 in Fig. 9 has the 
coordination number of 3 in a standalone square 16-point QAM constellation; but its 
coordination number is equal to 4 as a part of the 32-point QAM constellation. The 
ratio of this correction term to the value given by Equation (4.8) is of the order of the 
ratio of the perimeter of a 2 bj 12 x 2 bj 12 square to its area: 



correction term ~ 2 bj 12 1 2 bj -2 bj/2 —» 0, as bj —> oo 



Similarly, under the nearest-neighbor-error supposition, the first probability of the 
product in the second sum on the right hand side of Equation (7.1) is of the order of 

the ratio of the perimeter of a 2 bj 12 x 2 bj 11 square to the total number of constellation 
points 2 b , i.e., ~ 2 bj 12 1 2 b . The second - conditional - probability is, again, of the 
order of the ratio of the perimeter of a 2 bj/2 x 2 bj 12 square to its area, i.e., ~2~ bj!2 . 
Evaluating the whole expression in Equation (7.1) at large bj,j=l,..,J, yields: 



15 



J 2 bj 1 



1 + 0 



>b,/2 



f 1 > 



+ o 



(7.3) 



J 



20 



Using inequalities 1 < k < b t 1 2 and bj <b , j=l,..,J, and applying Equation (7.2) to 
Equation (7.3) yields: 



1 J 2 bj 1 J 2 bj 

p {k,b l ,..,b J ) = —T—+—Y— 



2 k %2 b 



2 k j^2" 



o 



v2 V2 



1 \\o 

2 k 



ryb 



1 + 0 



rtb/2 

J 



(7.4) 
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Finally, combining Equations (4.7) and (7.4) yields: 



p(k,b lt ..,bj) = -lj- 



f 1 \ 



1 + 0 



v2 J 



= P (k,bii+o 



( i 



\2< J 



(7.5) 



5 Details of the asymptotic structure of the p(k bj,..,bj) "tail" 

at k>bJ2 can be described in the following extension of Proposition 3. 

Proposition 4 . Let a non-square 2 h -point QAM constellation consist of 7/ 
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2 bl -point, J2 2 bl -point, and Ju 2 hu -point adjacent square sub-constellations, 



M 
m=\ 



and 



15 



b,<b 2 <..<b M <b 



Let each sub-constellation have been generated recursively according to the rule 
shown in Fig. 6. At large b m , m=l,..,M, 



20 



p{Kb x ,..,b J ) = p(k,b)\ + 0 



( 1 > 



2 b,2 J 



, as\<k<bJ2 



(7.6) 



p{k, b x ,. ., bj ) = p(k 9 b\ 1 - J x — 



v J 



f 1 A 



1 + 0 



J 



, asb 1 /2<k<b 2 /2 (7.7) 
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p(k,b l> ..,b J ) = p(k,b] I- J, 



— .- J. 



2 K\ 



2" m 2 b 



A 



1 + 0 



\2 j 



(7.8) 



asb m /2<k<b m+1 /2 



p{k,b x ,..,bj) = p(k,b)-J j 



2»m 



m 2 „ 



1 + 0 



( 1 ^ 



, asb M _J2<k<b M l2 (7.9) 



p(k,b l ,..,b J ) = 0 



v 2 y 



as k >b M 12 



(7.10) 
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Proof. 



Equation (7.6) has been proven in Proposition 3. 



At b m !2<k< b m+l 1 2 the first summation on the right hand side of 
Equation (7.1) is performed only over the sub-constellations with bit sizes larger than 

15 b m , that is ? in sub-constellations with smaller bit-size there cannot be two neighbor 
points with the Hamming distance larger than bJ2. At k > b M 12 p(k,b h ..,bj) in 
Equation (7. 1) is represented by the second sum only. There cannot be two neighbor 
points in the same square sub-constellation with the Hamming distance larger than 
by 2. Therefore, Equations (7.7) to (7.10) follow directly from Equations (7.3) and 

20 (7.4). 



Propositions 3 and 4 form a basis for approximating the main error 
correction properties (i.e., any probability distribution moment of interest) of a 
general, not necessarily square, constellation encoder by equations developed above 
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for square constellation encoders. This approximation can be used in error control 
calculations for QAM channels where both square and non-square constellation 
encoders are applied, for instance, in the G.lite/G.dmt ADSL standards. Since the 
series for determining moments (e.g., Equations (5.1) or (6.1)) are convergent, the 
5 contribution of the p(k f b If .. f bj) "tail" at k>bJ2 to these moments becomes 

negligible, as the sub-constellations sizes increase. The larger the constellation size, 
the more precise this approximation will be. 



8. Selecting a Bit Load 

10 

Referring to Fig. 10, a general flowchart of the present invention is 
shown. The modem initialization module 48 (Fig. 2) implements the method of the 
present invention. In step 182, a bit load is selected based on a self-similarity 
property of a constellation. In addition, one or more FEC parameters are also selected 
15 based on the self-similarity property. In step 1 84, data is transferred based on the bit 

load and FEC parameters. Alternately, the channel initialization module 43 (Fig. 2) 
implements the method of Fig. 10. 



In Fig. 1 1, a flowchart depicts an alternate embodiment of the 
20 invention that adjusts the target bit error rate to accommodate for fluctuations in the 

bit error rate. In step 192, a bit load is selected based on a self-similarity property of 
a constellation and by accommodating for fluctuations in the bit error rate. In 
addition, one or more FEC parameters are also selected. In step 194, data is 
transferred based on the bit load and FEC parameters. In one embodiment, the 
25 modem initialization module 48 (Fig. 2) implements the method of Fig. 1 1 ; 

alternately, the channel initialization module 43 (Fig. 2) implements the method of 
Fig. 11. 
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Fig. 12 depicts a more detailed flowchart of the present invention. 
Before describing the flowchart of Fig. 12, the relationships for determining the bit 
load and selecting FEC parameters will be described. 

Data transmission between two DMT modems uses a 
Reed-Solomon (RS) code for FEC in combination with a cyclic redundancy 
check (CRC) to detect and correct, if possible, errors. Data is transmitted in 
information frames. Each information frame of length N has an information field of a 
length of KRS symbols, a CRC field of a length of C symbols, and a redundancy 
field of a length of R symbols. The length of the redundancy field R determines a 
maximum number of symbol errors / that the Reed-Solomon code can correct. An 
information frame with t or fewer symbol errors will be corrected by FEC. An 
information frame with more than t symbol errors will pass through FEC with 
remaining errors. The CRC will then detect those errors and the frame will be 
retransmitted. 

Considering the DMT symbol as a superposition of n QAM 
waveforms, the size of a DMT symbol in bits B DM t^ equal to the sum of the number 
of bits carried by all useable sub-channels b t and is determined in accordance with the 
following relationship: 

i=\ 

The expression for the QAM symbol error rate ?j pg^, is determined as follows: 



^fe,yJ=i-t-(i-2^ /2 W4/ 3 - 1 ° r,/10/ l2 lk+1 -2)] 



\2 



(9.2) 
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where y t is the signal-to-noise ratio of the QAM waveform for the z th sub-channel. 
The uncoded bit error rate is determined in accordance with the following 
relationship: 



= a>[bAl-2 



~bJ2 



'■ ]erfc[j3-W- no /(2 b > +l -2)l2 - (l - 2 -*< y 2 )erfc{j3-W- no l{2 b >* -2)j 

(9.3) 

To determine the optimum bit load, Equation (9.4) is solved 
numerically for values of / and K. 



10 



\-W(t 9 K)sp = CD{b t (t 9 K% - T h ^ n )erfc{j3 . 10* 710 l{l b ^ 1 - 2)) 
J 

x [ 2 _ (i _ 2^ K)/2 )erf^ 



(9.4) 



where, W(t,K) - 



f K+C+R-l\ 
\ * ) 



-V{t + \) 



(9.5) 



15 In one embodiment, the value of co(b l (V, K)) is determined in accordance with 

Equation (5.4) of the present invention. Using Equation (5.4), the average fraction of 
erroneous bits per erroneous symbol (t, K)) is equal to 6){b i ) . 



20 



In another embodiment, in a system with automatic repeat request 
(ARQ) as well as FEC, the optimum bit load per sub-channel is determined in 
accordance with the following equation: 
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1- 



1 - W{t,K,k.)et l)h > j = coib^K,^ - 2 - b ^ KM)/2 )effc^3 -lO^ 10 /( 2 b ^ K ^ -2)) 
2-(l- T hXt ' K ^ n )erfc^3-W- no /{2 b ^ K ^ -2)1 



(9.6) 



where, W{t,K,ktr)- 



f K+C+R-l\ 

v * y 



1 _ _ ktr-l 

( K+C+R\~ 



(t+\)kt, 



V ^+1 J 



(t+\)hr 



(9.7) 



and k tr is equal to an allowable number of transmissions. Using Equation (5.4), 
fi)(& f (t, K, ktr)) is equal to a){b i ) . 

The SER level (RS symbol error rate) s s is related to a specified BER 
10 level s for the data at the output of the receiver after demodulation and applying error 
control in accordance with the following relationship: 



s s -1- 



(9.8) 



15 
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where a is the number of bits per RS symbol, and the coefficient J3 accounts for the 
effect of the descrambler. For a scrambler polynomial with m scr nonzero coefficients, 
the coefficient jB is approximately equal to the number of nonzero coefficients m scr 
ifi ~ wiser). Typically, the number of bits per RS symbol a is equal to eight; the target 
BER level sis equal to 10" 7 ; and the coefficient fi is equal to three. 

In another embodiment, to accommodate for fluctuations in the bit 
error rate, the target BER level s is modified by applying the mean square deviation 
using Equation (6.6) in which the target bit error rate sis equal to BER ' and the 
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adjusted target bit error rate s a dj is equal to BER 0 . The adjusted target bit error 
rate e a dj is used in Equation 9.8, rather than s , to determine the value of s s . 
Therefore, when solving Equations (9.4) and (9.6), the symbol error rate s s is 
modified to accommodate for fluctuations in the bit error rate. 

The net coding gain G n in bits/DMT symbol (the excess information 
transferred in one DMT symbol due to FEC) is determined in accordance with the 
following relationship: 



10 G n (t,K)^—^B DMT (t,K)-^-B DMT (0,K) (9.9) 



Equation (9.9) takes the FEC redundancy into account. The G.992.2 standard calls 
for the channel data rate to be transferred. Therefore the line coding gain G/ 5 that is, 
the total increase in the number of bits to be sent, is determined in accordance with 
1 5 the following relationship : 

^^^-^^^-^(O,^) (9.10) 

Equations (9.9) and (9.10) are modified when ARQ is applied. In another 
20 embodiment, when ARQ is applied, B DMT (t, K, ktr) is used rather than B DMT (t, K) . The 
net coding gain of a DMT symbol G n (t,K,k t r), that is, the excess information throughput 
per DMT symbol due to FEC and ARQ, is determined in accordance with the following 
relationship: 



25 



i *™ At ' KM) -^B DMT {W) (9.11) 



K+C+R v K+C 
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The average number of transmissions vis determined in accordance with the 
following relationship: 



v = 



1- 



K + C + R 
t + l 



1 


1- 


K+C+R 1 


lfktr 


s. 




J 




I t + l J 





(9.12) 



where k tr represents the maximum number of transmissions. 



10 



In yet another embodiment, since the procedure described yields 
fractional values for b t , an additional adjustment is applied. Such adjustment replaces 
the obtained b t values by their integer parts [p. J . Then a re-distribution of excess in 
bi over all sub-channels is performed. 



15 



A "mean-field approximation" for the number of bits of a DMT 
symbol in a channel is applied. The mean-field approximation is in accordance with 
Equation (9.13) as follows: 



n > v 



(9.13) 



i=l 



The effective number of sub-channels is represented by n e /f >? the effective 
20 signal-to-noise ratio for the channel is represented by y e ff, and the effective number of 
bits per sub-channel is represented by b(y e ff). The mean-field approximation reduces 
the complexity of the channel performance evaluation. Equations (9.4) and (9.6) 
perform equally well for the effective signal-to-noise ratio y eff as well as for the 
individual signal-to-noise ratios y of the subchannels. In Equation (9.13), the 
25 effective number of sub-channels n e f/ r and the effective signal-to-noise ratio #^are 

respectively determined in accordance with the following relationships: 
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n 



, and 



(9.14) 



n eff r,>r- 



(9.15) 



where y* is the threshold signal-to-noise ratio below which no information is passed 
through the channel Therefore, if the signal-to-noise ratio of a sub-channel does not 
exceed a predetermined threshold value y* for the channel, that sub-channel is not 
used. 



determination of the bit load and FEC parameters using the present invention and 
approximations (9.13) - (9.15). In one embodiment, step 202 of the flowchart of 
Fig. 12 is implemented in the modem initialization module 48 of the modem driver 40 
(Fig. 2); and steps 204-216 are implemented in the channel initialization module 43 
(Fig. 2). In step 202, one or more tables having a number of information bits to be 
supported by a sub-channel based on a self-similarity property of a constellation is 
generated for allowable sets of FEC parameters, and predetermined signal-to-noise 
ratios using the adjusted target bit error rate described above to accommodate for 
fluctuations in the bit error rate. In an alternate embodiment, the target bit error rate 
is used without adjustment. For each predetermined signal-to-noise ratio value y 9 
Equation (9.4) is solved at different allowable values of the maximum number of 
correctable errors t and the size of the information field K, and a number of bits per 
sub-channel; this solution, combined with Equation (9.13), yields the maximum 
number of information bits able to be supported by a DMT symbol within the 
system's constraints. In particular, when solving Equation (9.4), the average fraction 



Referring to Fig. 12, a more detailed flowchart depicts the 
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of erroneous bits per erroneous QAM symbol co(bi) is determined using 
Equation (5.4) which is based on the self-similarity property of the Hamming 
distances of a constellation. The tables are stored in memory. In an exemplary table 
shown in Table 2, below, a number of information bits to be supported by a 
5 sub-channel of a DMT symbol is shown for a set of parameters. 

In the channel initialization module 43 (Fig. 2), in step 204, 
measurements of the channel signal-to-noise ratio are made at different carrier 
frequencies associated with each sub-channel; and at the same time, a total number of 
10 '"working" sub-channels, that is, the effective number of sub-channels n eff , 

(Equation (9.14)) is determined. To measure the signal-to-noise ratio, the modem 
sends a predetermined signal pattern to its peer modem. In one embodiment, the 
predetermined signal pattern is a pseudorandom signal having a predefined spectrum. 

15 In step 206, an effective signal-to-noise ratio for the channel is 

determined based on the measured signal-to-noise ratio and the effective number of 
sub-channels n e /f using Equation (9.15). 

In step 208, a maximum net coding gain over allowable sets of FEC 
20 parameters is determined. In one embodiment, the tables generated in step 202 are 

iteratively accessed for each allowable set of FEC parameters based on the value of 
the effective signal-to-noise ratio to retrieve the associated number of bits per 
sub-channel. 

25 The net coding gain per sub-channel for each allowable set of FEC 

parameters is determined in accordance with Equation (9.16) below: 
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b(y eff ^K) 



K 



(9.16) 



n ejf K + C + R 



K + C 



In step 210, one or more FEC parameters are selected based on the 



maximum net coding gain per sub-channel g„. The selected FEC parameters yield the 
largest, i.e., the maximum, value of the net coding gain as determined in accordance 
with Equation (9.16). 



on the forward error correction parameters that provide the maximum net coding 
gain g n . The line coding gain per sub-channel gi is determined in accordance with the 
following relationship: 



In step 214, the bit load per sub-channel for the sub-channels, the line coding gain, 
and the one or more FEC selected parameters are transmitted to the peer modem. In 
step 216, data is transmitted using the selected bit load, line coding gain and FEC 
parameter(s). 



sub-channel b by SNR ^in decibels, and the number of correctable symbol errors t 
when s^lO" 7 , (3=3, a=8, C=0 and K=16 that is stored in memory. The values of the 
number of bits per sub-channel b are not the results of applying the average fraction 
of erroneous bits per erroneous symbol of the present invention, but are shown for 
illustrative purposes. 



In step 212, the corresponding line coding gain gi is determined based 




(9.17) 



In Table 2, below, an exemplary table depicts a number of bits per 
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Table 2: Exemplary portion of a table depicting a number of bits per sub-channel by 

SNR y in decibels, and the number of correctable symbol errors t when s=10~ 7 , 
P=3, q= 8, C=0, K=16 



y, dB 


b 


t=0 




10 


1. 023323e+000 


15 


2.078584e+000 


20 


3. 486503e+000 


25 


5. 066853e+000 


30 


6.711104e+000 


35 


8.376005e+000 


40 


1. 004683e+001 


45 


1. 171884e+001 


50 


1.339061e+001 


t=l 




10 


1.442432e+000 


15 


2. 674824e+000 


20 


4.183037e+000 


25 


5.808428e+000 


30 


7.473964e+000 


35 


9.151482e+000 


40 


l,083170e+001 


45 


l,251171e+001 


50 


1. 419067e+001 


t=2 




10 


1. 700372e+000 


15 


3.014947e+000 


20 


4 . 566292e+000 


25 


6.212193e+000 


30 


7.889203e+000 


35 


9.574868e+000 


40 


1. 126188e+001 


45 


1.294799e+001 


50 


1.463259e+001 
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Fig. 13 depicts an alternate embodiment of the channel initialization 
module 43 (Fig. 2), In this embodiment, the channel initialization module generates 
the bit load information, rather than the modem initialization module. Steps 204, 206, 
208-216 are the same as in Fig. 12 and will not be described. After performing 
5 step 206, step 218 determines a maximum number of information bits to be supported 

by a sub-channel for the effective signal-to-noise ratio, allowable FEC parameters and 
an adjusted target bit error rate based on a self-similarity property of a constellation. 
In an alternate embodiment, the target bit error rate is not adjusted. The bit load 
information is stored in memory in one or more tables similar to Table 2. 

10 

In another alternate aspect of the invention, these results can be 
applied to the optimization of DSL systems compliant with the G.992.2 (G.lite), and 
alternately the G.992.1 (G.dmt), standard. 



15 Although the present invention has been described with respect to a 

DMT communication system with multiple sub-channels, the present invention can 
also be used in a communication system having a single QAM channel 

Although various embodiments, each of which incorporates the 
20 teachings of the present invention, have been shown and described in detail herein, 

those skilled in the art can readily devise many other embodiments that still utilize 
these teachings. 



